Interaction induced localization in a gas of pyramidal molecules 
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We propose a model to describe a gas of pyramidal molecules interacting via dipole-dipole in- 
teractions. The interaction modifies the tunneling properties between the classical equilibrium 
configurations of the single molecule and, for sufficiently high pressure, the molecules become local- 
ized in these classical configurations. We explain quantitatively, without free parameters, the shift 
to zero-frequency of the inversion line observed upon increase of the pressure in a gas of ammo- 
nia or deuterated ammonia. For sufficiently high pressures, our model suggests the existence of a 
superselection rule for states of different chirality in substituted derivatives. 
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The behavior of gases of pyramidal molecules, i.e. 
molecules of the kind XY^ like ammonia NH^, has been 
the object of investigations since the early developments 
of quantum mechanics M. However, the behavior of these 
systems is still debated and some experimental facts re- 
main essentially unexplained B. 

For the single pyramidal molecule, owing to the great 
differences of characteristic energies and times, some 
adiabatic approximations hold. In particular, the one- 
dimensional inversion motion of the nucleus X across 
the plane containing the three nuclei Y can be separated 
from the rotational and vibrational nuclear degrees of 
freedom. The form of the effective potential for this mo- 
tion is a double well which is symmetric with respect to 
the inversion plane H. Due to tunneling across the fi- 
nite potential barrier, the eigenstates are delocalized in 
the two minima of the potential and, for energies below 
the barrier height, are grouped in doublets, i.e. couples of 
states with a relative splitting in energy small in compar- 
ison with the distance from the rest of the spectrum. For 
the pyramidal molecules under consideration, the ther- 
mal energy ksT at room temperature is much smaller 
than the distance between the first and the second dou- 
blet so that the problem can be reduced to the study of 
a two-level system corresponding to the symmetric and 
anti-symmetric states of the first doublet. 

The existence of delocalized stationary states is clearly 
in disagreement with the usual chemical view which, re- 
lying upon classical theory, considers the molecules as 
objects with a well defined spatial structure. In partic- 
ular, for the molecules under consideration the classical 
view predicts one of the two pyramidal configurations 
corresponding to the nucleus X localized in one of the 
wells of the inversion potential. 

The quantum prediction of stationary delocalized 
states implies the presence of a line in the absorption 
spectrum, the so called inversion line, at a frequency 
u — AE/h, where AE is the energy splitting of the first 



doublet. Experiments performed with NH3 Q, ND3 ||] 
and NT3 M reveal the existence of this inversion line 
in various rotational and vibrational bands. The fre- 
quency P of the inversion line has been measured as a 
function of the gas pressure P for NH3 [^ ^ and ND3 
Q. Starting from the expected value AE/h at P ~ 0, 
i^iP) decreases by increasing P and vanishes at a critical 
pressure Per — 1-7 atm for NH3 and Per — 0.1 atm for 
ND3. No quantitative theory has been proposed so far 
for this phenomenon. 

As early as 1949, in a short qualitative paper Q An- 
derson made the hypothesis that dipole-dipole interac- 
tion may induce a localization of the molecular states. In 
this way the important idea was introduced that inter- 
molecular interactions may be responsible for the ob- 
served phenomena. By considering a system of two or 
three interacting molecules, Margenau was able to pre- 
dict a decrease of the inversion frequency on reducing the 
distance between the moleculesjll . A quantitative dis- 



cussion of the collective effects induced by coupling a sin- 
gle molecule to the environment constituted by the other 
molecules of the gas was made in [|l2|. In this work it 
was shown that, due to the instability of tunneling under 
weak perturbations, the order of magnitude of the molec- 
ular dipole-dipole interaction may account for localized 
ground states. This suggests that a kind of phase transi- 
tion may be invoked to explain the behavior of NH3 and 
ND3 under variation of pressure. 

We have implemented this idea by constructing a sim- 
plified model of a gas of pyramidal molecules which ex- 
hibits the desired properties and allows a direct compar- 
ison with experimental data. Our model predicts, for 
sufficiently high inter-molecular interactions, the pres- 
ence of two degenerate ground states corresponding to 
the different localizations of the molecules. This transi- 
tion to localized states gives a reasonable explanation of 
the experimental results 0, |[ ^. In particular, it de- 
scribes quantitatively, without free parameters, the shift 



to zero-frequency of the inversion line of NH^ and ND3 
on increasing the pressure. 

We model the gas as a system of molecules nearly in- 
dependent in the following sense: each molecule is sub- 
jected to an external field representing the interaction 
with the rest of the gas to be determined self-consistently. 
We then analyze the linear response of this model to an 
electromagnetic perturbation to obtain the low energy 
excitation spectrum and its dependence on the inter- 
molecular interaction. Finally, we compare our theoreti- 
cal results with the available experimental data. 

We mimic the inversion degree of freedom of an isolated 
molecule with the Hamiltonian 



ho = -—a 



(1) 



where a^ is the Pauli matrix in the standard representa- 
tion with delocalized tunneling eigenstates 
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Since the rotational degrees of freedom of the single pyra- 
midal molecule are faster than the inversion ones, on the 
time scales of the inversion dynamics the molecules feel 
an effective attraction arising from the angle averaging 
of the dipole-dipole interaction at the temperature of the 
experiment ||I^. In the representation chosen for the 
Pauli matrices, the localizing effect of the dipole-dipole 
interaction between two molecules i and j can be repre- 
sented by an interaction term of the form crfcr^, where 
a^ has localized eigenstates 
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(3) 



In a mean-field approximation we obtain the total Hamil- 
tonian 



h{X) ^ -^a^ 
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(4) 



where |A) is the single- molecule state to be determined 
self-consistently by solving the nonlinear eigenvalue prob- 
lem associated to (0). The parameter G represents the 
dipolc interaction energy of a single molecule with the 
rest of the gas. This must be identified with a sum over all 
possible molecular distances and all possible dipole ori- 
entations calculated with the Boltzmann factor at tem- 
perature T. If g is the density of the gas, we have 
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where e^ is the relative dielectric constant, d the molec- 
ular collision diameter and the fraction in the inte- 
grand represents the Keesom energy between two clas- 
sical dipoles of moment ^ at distance r p3|. Equation 



(g) is valid in the high temperature limit which is ap- 
propriate for room temperature experiments. The exact 
expression, we assume g constant, is 



G = 
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where To = ji^ / {A-KeoSrd^kB) and In{u) is the modified 
Bessel function of the first kind. Expression (0) can be 
evaluated numerically and the result differs very little 
from (0) even for T <C Tg. For densities not too high, 
we set g — P/ksT so that, at fixed temperature, the 
mean-field interaction constant G increases linearly with 
the gas pressure P. By evaluating (^ we have 
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The solution of the eigenvalue problem associated to 
the Hamiltonian (Hh gives the following results. If G < 
Ai5/2, there is only one ground state |Ao) = |1), with 
energy 



En = —- 



AE 



(8) 



If G > AE/2, there are two degenerate ground states 



|A^) 


= V2+4gI^> + V2-4gI^> 


(9) 


|A?) 


= ^1A^), 


(10) 


with energy 
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By defining the critical value Gcr = Ai?/2, we distin- 
guish the following two cases. For G < Gcr, the ground 
state of the system is approximated by a product of de- 
localized symmetric single-molecule states corresponding 
to the ground state of an isolated molecule. For G > Gcr, 
we have two different product states which approximate 
the ground state of the system. The corresponding single- 
molecule states transform one into the other under the 
action of the inversion operator a'^ , see Eq. (p^), and, 
for G ^ Gcr, they become localized 



lim |A^) = \L) 

AE/G->0 
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The above results imply a bifurcation of the ground 
state at a critical interaction G = Gcr- According to Eq. 
(R) and using g = P/ksT, this transition can be obtained 
by increasing the gas pressure above the critical value 
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where Po = AE/d^. 

When the gas is exposed to an electro-magnetic radia- 
tion of angular frequency loq , we add to the Hamiltonian 
(m the perturbation 



hem{t) - efitW" 



(14) 



where e is a small parameter and f{t) =^ 9{t) cos{LUot), 
6{t) being the Heaviside function. The choice of a dipole 
coupling approximation, h^m oc cr^, is justified for a ra- 
diation of wavelength long with respect to the molec- 
ular size. Under the effect of the perturbation ( 14 ) the 
single-molecule state |A(f)) evolves according to the time- 
dependent nonlinear Schrodinger equation 
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^[h(X(t)) + ef{t)a^]\X{t)), (15) 



with h{X) given by (H). The linear response to the pertur- 
bation (|^ is expressed by the generalized susceptibility 
Q TZ{lu) = Si{uj)/f{uj), where f{uj) and Si{uj) are the 
Fourier transforms of f{t) and Si{t), with Si{t) defined 
by 



S{t) EE {X{t) I a^ I X{t)) 
= So{t) + e Si{t) + 



(16) 



Let us assume that at time i = each molecule is in 
the delocalized ground state |Ao) = |1). The solution of 
Eq. ( p^ ) with the initial condition |A(0)) = |1), gives 



n{Lu) = 



2AE 



[hujY ~ {AE^ - 2GAE) ■ 



(17) 



The generalized susceptibility has a unique pole at pos- 
itive frequency which corresponds to the inversion line 
frequency 
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The residue of TI{lo) at this pole, namely 
(1 — 2G/AE) ' , represents the corresponding transi- 
tion probability. 

Now we compare our theoretical analysis of the inver- 
sion line with the spectroscopic data available for am- 
monia ItI pi and deuterated ammonia 0. In these ex- 
periments the absorption coefhcient of a cell containing 
NH^ or ND^ gas at room temperature was measured at 
different pressures. The resulting data are reported in 
Fig. n^. The frequency v of the inversion line decreases 
by increasing P and vanishes for pressures greater than 
a critical value. There is factor about 15 between the 
critical pressures of NH3 and ND3. 

By using Eq. (R) with g = P/ksT, the theoretical 
expression ( [l8|) for the inversion line frequency becomes 
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(19) 
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FIG. 1: Measured inversion-line frequency ;? as a function of 
the gas pressure P for NHz (dots, left and bottom scales, 
data from 10, M]) and ND-^ (squares, right and top scales, 
data from B]). The solid line is the theoretical formula (119) 
with Per = 1.695 atm for NHi and Per = 0.111 atm for NDi 
calculated according to (|l5 



where Per is given by (|l^). Note that this expression 
does not contain free parameters. We used the follow- 
ing values taken from |, |l|: ^ = 1-47 D, d = 4.32 A, 
APjvifg = 0.81 cm"\ AEnd-s = 0.053 cm"^ Assuming 
£r — \ and T = 300 K, we obtain Per = 1.695 atm for 
NH^ and Per = 0.111 atm for ND^. The agreement of 
the theoretical v{P), also shown in Fig. 1, with the ex- 
perimental data is impressive considering the simplicity 
of the model. 

Equation ( p9| ) predicts that, up to a pressure rescal- 
ing, the same behavior of v{P) is obtained for different 
pyramidal molecules 
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where 



= Pe, 



/Pc, 



In the case of ND'i 



7 — ^crX'Yg'/ ^ciXYs- 

and NH^, at the same temperature T we have 7 = 
AEmh^/ AEnds — 15.28. This factor has been used to 
fix the scales of Fig. 1. We see that in this way the NH^ 
and ND^ data fall on the same curve. 

The intensity / of the inversion line predicted by our 
theoretical analysis is given, up to a constant, by the 
product of the photon energy hv and the residue of (|l^). 
The divergence of the transition probability is cancelled 
by the vanishing photon energy and we obtain for a gas 
of N molecules / oc NAE. By writing N = PV/ksT, 
V being the volume of the absorption cell containing the 
gas, we have that at fixed temperature the line intensity 
increases linearly with pressure. This behavior is con- 
firmed by the experimental data 0, §, ^ • 

It is interesting at this point to compare our approach 
with a previous study [Rra where the experimental data 



are successfully reproduced by a formula with three free 
parameters for the shape factor of the inversion line. In 
pa these parameters are determined by trial and error. 
In our work we propose a simplified theory for the shift of 
the inversion line based on the dipole-dipole interaction 
and there are no free parameters. The specific predic- 
tion of our model for the critical pressure Per in terms of 
the electric dipole n of the molecule, its size d, the split- 
ting AE and the temperature T of the gas, successfully 
verified in the case of ammonia, should be experimentally 
tested also for other pyramidal gases. We emphasize that 
the study of the inversion spectra of pyramidal molecules, 
like NH3 and ND3, in recent years has acquired a con- 
siderable interest in geophysical |1^ and astrophysical 
research Q. 

It is not easy to assess the region of validity of the 
model starting from first principles, i.e. considering all 
the degrees of freedom of the molecules. It is reasonable 
to assume that the mean field calculation is meaningful 
as long as the interactions among the molecules are small 
compared to the inversion line frequency ([l8|). If we take 
the model seriously, at least for the qualitative aspects, 
at pressures greater than the critical one we have the 
following situation. In the limit of an infinite number 
of molecules, the Hilbert space separates into two sectors 
generated by the ground state vectors given in mean field 
approximation by 






|Ao^> . 



(21) 
(22) 



These sectors, which we call Hl and Hr, cannot be con- 
nected by any operator involving a finite number of de- 
grees of freedom (local operator). According to g] this 
means that a superselection rule operates between the 
two sectors distinguished by the eigenvalues of an observ- 
able. It is natural |19] to define the chirality operator 



1 ^ 
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It is immediately verified that 
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(24) 



for tjj in Hl and Tin, respectively. This mean value dis- 
tinguishes the two sectors but is a property weaker than 
■0 being an eigenstate of x- According to (n2|) only in 
the limit AE/G -^ 0, the states ip become completely lo- 
calized and therefore eigenstates of x- The fact that the 
ip are not eigenstates of x for AE/G 7^ is connected 
to the non-orthogonality of the mean-field one-molecule 
states jAg) and |A^) and is presumably an artifact of the 
approximation. 

This analysis has an interesting implication. Our 
model applies not only to molecules XY^ but also to their 



substituted derivatives XYWZ. In fact, the difference 
should not matter as far as the shift of the inversion line 
is concerned. However, an important difference between 
the two cases is that for XY^ the localized states can be 
obtained one from the other either by rotation or by space 
inversion, while for XYWZ they can be connected only 
by space inversion. This implies that XYWZ molecules 
at a pressure greater than the critical value are chiral 
and therefore optically active. Recent ultrasensitive ex- 
perimental methods have been developed which allow the 
quantitative measurement of optical rotation in gaseous 
compounds |2^. It would be interesting to measure the 
optical activity of pyramidal gases for P > Per- 
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